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Abstract 

Using high-precision numerical analysis, we show that 3 + 1 dimensional gauge the- 
ories holographically dual to 4 + 1 dimensional Einstein-Maxwell-Chern-Simons theory 
undergo a quantum phase transition in the presence of a finite charge density and 
magnetic field. The quantum critical theory has dynamical scaling exponent z = 3, 
and is reached by tuning a relevant operator of scaling dimension 2. For magnetic field 
B above the critical value the system behaves as a Fermi liquid. As the magnetic 
field approaches Be from the high field side, the specific heat coefficient diverges as 
1/{B — Be), and non- Fermi liquid behavior sets in. For B < Be the entropy density 
s becomes non- vanishing at zero temperature, and scales according to s ^ y/Bc — B. 
At B = Be, and for small non-zero temperature T, a new scaling law sets in for which 
s ^ T^/^. Throughout a small region surrounding the quantum critical point, the 
ratio s/T^/^ is given by a universal scaling function which depends only on the ratio 

(S-Se)/r2/3. 

The quantum phase transition involves non-analytic behavior of the specific heat 
and magnetization but no change of symmetry. Above the critical field, our numer- 
ical results are consistent with those predicted by the Hertz/Millis theory applied to 
metamagnetic quantum phase transitions, which also describe non- analytic changes in 
magnetization without change of symmetry. Such transitions have been the subject of 
much experimental investigation recently, especially in the compound Sr3Ru207, and 
we comment on the connections. 



^This work was supported in part by NSF grant PHY-07-57702. 



1 Introduction and summary of results 



The AdS/CFT correspondence provides a precise and powerful tool for the study of ther- 
modynamics, statistical mechanics, and transport properties in a variety of 4-dimensional 
gauge theories at finite temperature, charge density, and magnetic field. In the large and 
large 't Hooft coupling limits these gauge theories are holographically dual to certain black 
brane solutions (with both electric and magnetic charges) to 5-dimensional Einstein-Maxwell 
theory. This theory includes a Chern-Simons term, whose coupling k captures the strength 
of the chiral anomaly in the Maxwell current. For the special value k = 2/^/3 (in our 
conventions) the bulk Einstein-Maxwell theory is a consistent supersymmetric truncation of 
Type IIB or M-theory [H El |3j, while for k not equal to this value, supersymmetry is lost. 

The study of thermodynamics and transport properties in 4-dimensional Yang-Mills the- 
ory is of direct interest to the physics of heavy ion collisions at RHIC (and soon at the 
LHC), where quarks and gluons are subject to high temperatures, strong magnetic fields, 
and large currents. Of great interest as well is the possibility of studying novel phases of 
finite density matter at low temperatures. By varying the charge density p and magnetic 
field B one can search for signals of zero temperature quantum phase transitions jlj, as is 
done experimentally in, for example, the heavy fermion compounds and high temperature 
superconductors . 

The thermodynamics of the relevant solutions to 5-dimensional Einstein-Maxwell theory 
can be understood analytically in two limiting cases. In the absence of magnetic fields, the 
relevant supergravity solution is the electrically charged (Reissner-Nordstrom) black brane in 
AdS^. At high temperature, its entropy density s scales as , as dictated by scale invariance 
in the UV, and in accord with weak coupling Yang-Mills theory. As T tends to 0, the black 
brane tends to its extremal form with AdS2 x near-horizon geometry, and non- vanishing 
entropy density. On the one hand, this AdS2 factor is needed in "semi-holographic" models 
of non- Fermi liquid behavior O |6l [3 EJ (see also [9j). On the other hand the ground state 
entropy density is rather exotic from the Fermi surface perspective. It is exotic also from 
the point of view of the CFTs that arise in the AdS/CFT duality, which typically contain 
massless charged bosons that would be expected to condense and lead to a unique ground 
state. 

The other limiting case consists of vanishing electric charge density and finite magnetic 
field B. At high temperature the entropy density s again scales as T^, but at low temperature 
the analysis of [10] showed that s ^ BT. The overall numerical coefficient in this linear 
scaling law was computed analytically by taking advantage of the near horizon AdSa x 
factor that emerges in this regime. On the CFT side, the low temperature physics was seen to 
be controlled by a gas of fermions arising from the lowest Landau level of the 4-dimensional 
gauge theory in the presence of a magnetic field. 

The general case of nonzero B and p was studied in [Hj, and the low temperature 
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thermodynamics were found to depend crucially on the value of the Chern-Simons coupling k. 
This was characterized in terms of flows in parameter space as the temperature was lowered. 
These flows head towards three distinct flxed points, depending on the value of k. For 
< 1, the geometry near the horizon can be thought of as AdS2 x deformed by the 
presence of the magnetic fleld. The entropy density was found to be non- vanishing at T = 0. 
Precisely at A; = 1, the solutions support a near-horizon warped AdSa x factor (the 
warped solutions were studied in the context of topologically massive gravity in [121 fT3]). 
and the entropy density is again non- vanishing at extremality. For k > 1 (which includes the 
supersymmetric value k = 2/^/3), the presence of a moderate strength (or larger) magnetic 
field was seen to lead to a precipitous drop in the entropy density as the temperature was 
lowered. However, our numerics were found to break down in the combined regions of very 
low temperatures and small magnetic fields, presumably as a result of one of our choices of 
gauge and, as a result, our analysis stopped short of obtaining reliable data for the entropy 
density at low T over the full range of magnetic fields. 



1.1 Summary of results 

In the present paper, we shall apply a simple remedy to the gauge choice problem which 
plagued the low temperature numerical work of , and carry out high-precision numerical 
analyses down to ultra-low temperatures, for a wide range of values of / p^. 

By virtue of scale invariance, only dimensionless combinations of quantities, such as 
I (? ^ afltord any intrinsic physical meaning. Thus, we shall introduce normalized, dimen- 
sionless, magnetic fleld 5, temperature T, and entropy density s via the following relations, 

R = A f- ^ ' nn 

- p2/3 - (^3 + p2)l/6 ^ ~ (53 + p2)l/2 V • ) 

The results are displayed schematically in Figure [T] and summarized below. 

1. Our main result is that a continuous quantum phase transition occurs at 

5^ = 0.4994240 ± 0.0000007 (1.2) 

2. For B > Bc^ the entropy density goes to zero linearly with temperature, s ^ 
reflecting the same Fermi liquid physics that was seen in the large B limit. Here, 
however, the linear behavior is occurring at flnite charge density where the result does 
not follow from conformal invariance. Note that the speciflc heat C, at constant 5, 
has the same behavior as the entropy density, since C = Tds/dT ^ T. 

3. On approaching Be from the large B side, the coeflScient of the linear term in the s 
versus T relation is found to diverge according to, 

^ ^ a = 1.003 ± 0.005 (1.3) 



T (B- B,Y 



3 




Figure 1: Schematic phase diagram iUustrating the various behaviors of the entropy density 
versus temperature and magnetic field. The region inside the dotted hne is controUed by 
the quantum critical point at (T = 0, 5 = 5c), and the entropy density can be expressed in 
terms of a sing: le scaling function / of (5 - B^jT'^l^. We move around inside this region by 
changing the temperature T and the relevant coupling B — Be- The boundary of the region 
is defined to be where irrelevant operators become important. The yellow region denotes a 
regime where temperature is the largest energy scale, corresponding to the argument of the 
scaling function / being small. Outside the yellow region the low temperature behavior of 
the entropy density, for fixed 5, is either constant or linear in T, depending on whether the 
quantum critical point is approached from below or from above as T ^ 0. 



This signals a breakdown of Fermi liquid behavior. 

4. Approaching the quantum critical point along the temperature axis at fixed B = Bc^ 
the entropy density exhibits a new power law scaling, 

a = 0.335 ± 0.005 (1.4) 

5. For B < Be the zero temperature entropy density is found to be nonzero. Near Be it 
obeys the scaling, 

{Be- BeY T = 0.500 ± 0.001 (1.5) 
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6. The entropy density in the vicinity of the fixed point can be expressed in terms of a 
scahng function / as, 




This is to be contrasted with regions far from the critical point, where the entropy 
density is a nontrivial function of two dimensionless combinations of p, B and T. 

From our results, we can infer that the quantum critical theory lives in 1 + 1 spacetime 
dimensions (that is, there are no long range correlations in the remaining 2 spatial directions 
along the boundary), has dynamical critical exponent z = 3 (since the value found numeri- 
cally for a is consistent with a = 1/3), and has a relevant operator with scaling dimension 2. 
The relevant operator corresponds to a change of B away from Be. 

1.2 Comparison with known quantum critical systems 

It is illuminating to place these results within the context of known quantum critical systems. 
Although thermodynamic quantities such as the specific heat and magnetization behave in a 
non-analytic fashion across the phase transition, we note that there is no change of symmetry 
associated with the transition. A finite temperature example of such behavior is the liquid- 
gas transition in water. A zero temperature version is a metamagnetic quantum critical 
point p!4i ITS] , whose behavior closely parallels our system (a finite temperature version 
of a holographic metamagnetic phase transition in the D4-D8 system was studied in [16\). 
Metamagnetism refers to a sharp change in the magnetization of a material as an external 
magnetic field is tuned through some nonzero value. At finite temperature, metamagnetic 
transitions typically consist of first order lines terminating at a second order critical point, as 
in the liquid- vapor case. If the critical point can be brought to zero temperature by adjusting 
some control parameter, one obtains a metamagnetic quantum critical point. 

The standard approach to such a quantum critical point is based on the Hertz/Millis 
theory [TTl [IHl [19], for which the effective action in momentum space is, 

S = |rfcjd^fc(^^ + A;2 + (S-4)^ |0(c^,fc)|2 + ... (1.7) 

The real bosonic field cj) represents the local magnetization, and the above action can be 
obtained by integrating out fermions at 1-loop. Under scale transformations acting as ^ 
A/c, we see that uj should be assigned scale dimension 3, and hence z = 3. Similarly, B — Be 
has scale dimensions 2. These assignments match what we found for our system, which 
furthermore corresponds to d = 1, since the Landau level quantization only allows low 
energy modes to propagate parallel to B. 
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The action (1.7) is only meant to be applied for B > Be. Indeed, the behavior of our 
system in the region B < B^^ with its nonzero ground state entropy density, cannot be 
described by this action alone. 

Metamagnetic quantum criticality in the compound Sr3Ru207 has been the subject of 
extensive experimental investigation in the past few years (see [15j and references therein), 
and we will comment on this connection in Section 3.8. 

A large number of other AdS/CFT examples undergoing phase transitions have been 
studied in the literature, both continuous and discontinuous, and at finite and zero temper- 
ature; for example [20l[2ll|22l[23l[24l|25ll^ However, our setup seems particularly 
attractive and is nicely related to real experimental systems. In particular, unlike other ex- 
amples of quantum phase transitions we do not have to add any extra ingredients in the way 
of scalar fields or probe branes. We employ only a metric and an Abelian gauge field, with an 
Einstein-Maxwell-Chern-Simons action that is known to describe all supersymmetric AdSs 
theories related by compactification of Type IIB or M-theory [il[2[|3]. Thus our framework 
is both simple and universal. 

The remainder of this paper is organized as follows. In Section 2, we spell out the set-up 
of the holographic calculations, including the specification of initial data at the horizon, 
asymptotic data at the AdS^ boundary, and the construction of regular gauge choices. In 
Section 3, we present the results of high-precision numerical solutions to the reduced Einstein- 
Maxwell-Chern-Simons equations, identify the quantum critical point, and critical behavior 
in its vicinity. We also compare our results with those from the eflFective Hertz/Millis theory, 
and comment on the connection with metamagnetic quantum criticality observed in real 
materials like Sr3Ru207. A discussion of the results and open avenues for future research is 
given in Section 4. 



2 Holographic calculations 

In this section we spell out some of the technical details of our computations. Results of the 
numerical calculations will be presented in section 3, and the impatient reader may wish to 
jump there. 

The starting point for our holographic calculations is 5-dimensional Einstein-Maxwell 
theory with a Chern-Simons term. Throughout this paper, the Chern-Simons coefficient 
k will be considered fixed at its supersymmetric value k = 2/^/3. A detailed discussion 
of the action, including boundary terms needed for holographic renormalization, and the 
construction of the boundary current and stress tensor may be found in pjj. Here, we shall 
limit our discussion to the field equations and the asymptotic behavior of the fields near the 
horizon, and at the asymptotic AdS^ boundary. The Einstein-Maxwell field equations are 
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given by dF = and, 



= d^F + kFAF 
Rmn ^Qmn + -F^^FpQ qmn "^FmpF^^ (2.1) 

Uniformity and constancy in time of the magnetic field B and the charge density p aUow us 
to restrict to a space-time translation invariant Ansatz, given by, 

F = Edr Adt + Bdxi A dx2 + Pdxs A dr 
ds^ = U-Hr^ -Ude + e^^ (dxl + dx^^+e^"^ {dx^ + Cdtf (2.2) 

The functions E^P^U^ C depend only on the radial coordinate r, while the magnetic 
field B is constant by the Bianchi identity. A gauge choice has been made here for the 
coordinate r in order to put the Ansatz in canonical form with matching coefficients of its 
first two terms in ds^. The reduced field equations were given in [llj. 

2.1 Data at the horizon 

The reduced field equations are to be solved subject to regularity conditions at the horizon 
and at the asymptotic AdS^ boundary r = oc. For the purpose of numerical analysis, it 
will be convenient to parametrize solutions in terms of data at the horizon which satisfy 
the regularity conditions at the horizon. Regularity of the full solution, including at the 
asymptotic AdS^ boundary, must then be verified numerically for each set of data. 

We begin by spelling out the data at the horizon. (This discussion will parallel the one 
presented in [11], but there will be important differences motivated by the need to remedy 
the gauge choice problems alluded to in the Introduction.) The (outer) horizon at r = r+, 
and the Hawking temperature T are defined by, 

[/(r+) = 47rr = U\r^) (2.3) 

Our numerical analysis will always be carried out at T 7^ 0, even though T may become 
very small; thus, we are free to rescale t, and set U'{rj^) = 1. By rescaling also Xi,X2,X3, 
we may set V{r^) = W{r^) = 0. Invariance of the Ansatz under asymmetry^ (under which 
Xs ^ Xs — at^ C ^ C + and E ^ E — aP)^ allows us to set C(r+) = 0. With these 
choices, the fields at the horizon take the form, 

Fh = qdr Adt + b dxi A dx2 + p dxs A dr 

ds]j = dxl + dxl + dxl (2.4) 

The reduced field equations of [llj for P, Vl^, C, combined with the requirement of 
regularity at the horizon, dictate certain relations amongst the remaining data at the horizon. 
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namely q,b,p,V'{r-^), VF'(r+), and C"(r+). They are given by (for U'{r+) = 1), 



= 
3V'(r+) = 



p-g(C"(r+)- 
12 - - 45^ 



6W^'(r+) = 24-4{q'^ -b'^) -3C'{r+f 



(2.5) 



Given (6, g), the data p and C(r+) are related to one another by the first equation of (2.5). 
Keeping either p or C(r+) as an independent free parameter, the remaining initial data 



V\r^) and W\r^) are uniquely determined by the last two equations of (2.5). 



2.2 Gauge fixing and regularity 

The parametrization of the horizon data presented above is more general than the one given 
in [15, since here the parameter p is kept unspecified, while p was set to in [11]. The 



argument invoked to set p = was the covariance of the Ansatz of (2.2) under boosts in 
the X3-direction. Under a boost by velocity /3, the space-time coordinates transform as 
usual, t ^ i = j{t — /3xs) and Xs ^ Xs = 7(^3 — /3t) with 7^(1 — 1, and where /3 

cannot exceed the speed of light, \/3\ < 1. The transformation under boosts of the functions 
{E^ P, [/, C) P, [/, y, (7), however, must be accompanied by a transformation 

of the holographic coordinate r ^ which is required to restore the boosted Ansatz back 



to the canonical form of (2.2). As a result, the Maxwell fields transform as follows, 

E{r)dr = -i(E{r) - pP{r))dr 

P{f)df = -i(P{r) - l5E{r))dr (2.6) 

where U{f)~^df = U{r)~^dr. The ratio p/q transforms diS p/q = {p/q — /3)/(l — f3p/q). 
Clearly, provided > p^, a boost hy /3 = p/q may be used to set p to zero. The problem, 
however, is that the coordinate transformation r ^ r required to accompany this boost may 
be singular on some of the functions E^ P^U^V^W^C. For example, the transformation law 
under a boost by velocity /3 of the function C is given by, 

G (r) ^ G (r) - (i + pdrWe^^i^) - P^U{r) ^^"^^ 
We have verified numerically that, in the region of low T where our numerics were found to 



break down in [llj, the denominator in (2.7) indeed crosses zero as r is increased away from 
the horizon at r+, thereby rendering this boost transformation singular. 

The remedy to this problem is simple: the parameter p should be left unspecified, thereby 
eliminating the need to perform boost transformations on the fields. It will be convenient 
to parametrize solutions by the values 6, g, C(r+), so that the remaining initial data at the 
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horizon, y(r+), and W\r^) are uniquely determined by (2.5). Not every assignment of 
b^q^C^{r^) will produce a regular solution. Also, two regular solutions may be related to 
one another by a regular boost, and thus be physically equivalent. The parameter space of 
all regular solutions may be described as follows. To every pair (6, q) we assign the maximal 
interval r(6, g) on the real line such that for every value of C\r^) G r(6, g), the solution 
specified by the parameters b^q^C\r^) is regular. The end points of the interval r(6, g) 
correspond to the limits where the velocity of the solution tends to the speed of light. When 
computing boost invariant physical quantities, such as the magnetic field S, and the rest 
frame temperature T, charge density p, and entropy density 5, the datum C(r+) may be 
chosen to be any value in the interval r(6, q). 



2.3 Data at the asymptotic AdS^ boundary r ^ oc 

For regular solutions the functions [/, V^W^E^P^C have the following asymptotic behavior 
as r ^ oc, (keeping only leading contributions), 

E ^ {es- coPsy-^ P - psr-^ C cq (2.8) 

In these coordinates, the conformal boundary metric is —dt'^ + v{dxl + dx^) +w{dxs + Cgrft)^, 
and the solution's velocity is ^/w Cq. Rescaling Xi, X2 by y^, and Xs by y^, while combining 
an ^-transformation with a boost by a velocity /3 = —^/wco^ restores the coordinates to the 
standard Minkowski metric, and yields the following expressions for the physical magnetic 
field S, temperature T, charge density p, and entropy density s in the rest frame of the 
soiufcioiajj (see [llj for derivations). 



B=^ T=]^ p = 7c(e3 - coPs) s = (2.9) 



Here, 7^ = (1 — wcl) and the normalized entropy density s is defined as 5 = GsS'/Vol, 
where S is the total entropy in volume Vol, and G5 is the 5-dimensional Newton constant. 

By virtue of scale invariance, only dimensionless combinations of quantities, such as 
B^/p^^ afford any intrinsic physical meaning. Thus, the results on phase transitions and 
associated critical exponents at critical points, aimed for in this paper, will all be derived 
by evaluating the dependence between the dimensionless physical quantities s and T, for 



various (fixed) values of 5, which were all defined in (1.1). 



^In the CFT, the rest frame corresponds to a statistical ensemble weighted by the Boltzmann factor 
^-{H-iiQ)/T ^ Boosting produces additional chemical potentials multiplying momenta and currents. 
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2.4 Numerical fine points 

The initial data for any regular solution is a pair (6, q) and a value C(r+) G r(6, g), where 
the interval r(6, g) was defined so that the data 6, C(r+) produce a regular solution. For 
certain values of (6, g), the interval r(6, q) may be empty. For k = 2/^/3, numerical analysis 
yields the following bound, 

r(6, g) = when + 26^ > 6 (2.10) 

The precise form of the critical curve C in the (6, g)-plane inside of which r(6, g) 7^ is not 
known analytically, but may be obtained numerically, as was done in [TT] . 

Away from the low temperature regime the interval r(6, q) is sufficiently large that it 
is possible, at least for k = 2/\/3, to make a single uniform choice for C^(r+) and still 
cover most of the parameter space. A convenient choice is C(r+) = 26. However, as the 
temperature is lowered r(6, q) shrinks, and one is forced to tune C^(r+) to greater precision. 
This effect becomes especially pronounced in the region of low magnetic field where a ground 
state entropy develops. However, by specifying more general tunings for C(r+) as we vary 
(6, q) we are able to fully cover the low temperature region. 

The approach to ultra-low temperatures, which is needed in various parts of our numerical 
work, requires a high degree of fine-tuning of the horizon initial data (6, q) and the gauge 
choice C(r+). It also requires evaluating the asymptotic data, such as Cq, 63 and ^3 at 
large values of r, which we typically have taken to range from 10^^ to 10^^. With such high 
degrees of fine-tuning, and extended ranges of integration, the issue of numerical accuracy 
and numerical stability of the calculations becomes of utmost importance. Our calculations 
were performed with 15 to 20 digits accuracy, and the ODEs were solved with absolute and 
relative error tolerances ranging from 10~^° to 10~^^ for the lowest temperatures. Stability of 
the results was checked versus changing the asymptotic value of r, the absolute and relative 
error tolerances, and the number of digits. 

3 Numerical results 

In this section, we present our numerical results, organized as a function of the magnitude 
of the magnetic field 5, starting at large B. 

3.1 Large B regime 

In [11] the case B ^ .53 (quoted there as B^ ^ .15) was studied at the supersymmetric value 
k = 2/\/3. As displayed in Figure 3 of pTj the low temperature entropy density was seen 
to drop well below its S = value, and appeared to be heading towards zero. However, the 
numerics broke down at temperature T .02, and so did not allow for an exploration of 
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ultra-low temperatures. By adjusting C(r+) as we lower the temperature, we can do much 
better, as shown in Figure [2] below. The entropy density clearly vanishes at zero temperature. 
A linear behavior, 5 ^ T is evident in the approach to zero temperature, as shown in the right 
panel of Figure^ By repeating this analysis for other sufEciently large values of S, we find 




Figure 2: Plot of entropy versus temperature. On the left we compare 5 = .53 to 5 = 0; 
this plot is an improved version of Fig. 3 in [11]. On the right we exhibit the linear s ^ T 
low temperature behavior. 

similar behavior: the entropy density vanishes at zero temperature, and does so linearly at 
very low temperatures. This linear behavior is characteristic of Fermi liquids. The textbook 
intuitive explanation for the linear dependence is that at low temperature the smoothed out 
step function form of the Fermi-Dirac distribution implies that only electrons within energy 
ksT of the Fermi energy contribute. We indeed expect our system to be described by a 
theory of light fermions in this large 5, low T regime, since the magnetic field raises up all 
energies except for the lowest fermion Landau level. Our numerical results are a pleasing 
confirmation of this intuition. 

In the B ^ 00 limit it was found in [10] that a near horizon AdS^ x factor emerged 
at low energies, and the resulting 1 + 1 dimensional conformal invariance could be used to 
explain the s ^ T behavior. This is no longer true away from this limit, as the charge density 
p introduces an additional scale into the problem, and from the numerics we can see that 
the near horizon geometry becomes deformed away from AdS^ x R^. 

Charged black hole solutions involving scalar fields whose entropy densities similarly go 
to zero at extremality have been studied in [29l [30l [3ll [32] . 



11 



3.2 Approaching Be 

Experimentally, a breakdown of Fermi liquid behavior upon tuning an external magnetic 
field can often be seen in a divergence of the specific heat coefficient, defined as, 
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C 
f 



(3.1) 



Since C = Tds/dT^ we can equivalently write 7 = s/T at low temperatures. In our system, 
7 stabilizes at a constant value for large S, but is seen to diverge at a critical value B = Be. 
Numerically, the critical value Be is found to be bounded as follows. 



0.124568 < Bt < 0.124569 



(3.2) 



which results in the the value quoted already in ([L2|), namely Be = 0.4994240 ± 0.0000007. 
To characterize the divergence we display a plot of s/T versus 1/{B — B^)^ as T ^ in 
Figure [3) The straight line shows that the low temperature entropy behaves as 



T B -Br 



(3.3) 



as T ^ for fixed B near (but larger than) Be. Upon approaching Be we find that the linear 
regime with 5 ^ T is confined to an ever smaller low temperature window. 



s_ 

f 




B-Br 



Figure 3: Plot showing the divergence of s/T near Be at low temperatures. The straight line 
through the data points is included to guide the eye. 
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3.3 Scaling at the critical magnetic field B = Be 



Next we set B = Be and again study the low temperature behavior of the entropy density. 
We find a new scahng law, s ^ T^/^, which numerically extends over at least four orders of 
magnitude in temperature, 10~^ < T < 6 x 10~^, as shown in Figures [4| This nontrivial 
power law manifestly represents non- Fermi liquid behavior, analogous to what is seen in real 
materials at a quantum critical point. 

Towards ultra-low temperatures, the numerical behavior of s will ultimately turn over to 
be linear in T (for B > Be) or to be a non-zero constant for B < B^ (see the next subsection). 
This deviation fro m s ^ T^/^ scaling is caused by the fact that the value of B^ is known only 
numerically from (1.2), and we are never able to sit at precisely the value of B. 



— In s 



B = Br 



hiT 



Figure 4: Plot showing s ^ T^/^ scaling behavior at 5 = Be. The straight line through the 



data points is included to guide the eye, and has slope 1/3, consistent with (1.4). 



3.4 The quantum critical region and crossover 

The s ^ T^/^ scaling law extends to the vicinity of the critical magnetic field 5c, but the 
range over which this scaling law holds shrinks as \B — B^ is increased. The corresponding 
behavior of s is illustrated in Figure 5. 

First, on the right panel of Figure 5, — ln{s) is plotted versus — ln(T) for B^ = 0.124569. 
The T^/^ scaling law exhibited in the preceding subsection, is clearly recovered here. At 
sufficiently high temperature, the T^/^ behavior eventually crosses over to the dependence 
controlled by the UV theory. This is clearly shown on the right panel of Figure 5, where the 
cross-over region may be identified with the temperature interval 0.02 < Tcross-over < 0.5. 
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Second, on the left panel of Figure 5, the flows of — ln(5) as a function of — ln(T) at 
various flxed values of B are shown. Curves a, b, c, d, e, and f clearly exhibit the turnover 
from s ^ T^/^ scaling behavior to linear s ^ T behavior at ultra-low temperatures. Given this 
cross-over behavior, we know that the corresponding values of B must all be above the critical 
magnetic field 5c, with value closest to critical corresponding to curve f with B^ = 0.124569. 
Curve g behaves completely differently at ultra-low temperatures, and s is seen to tend 
towards a non-zero constant. Given this behavior, we know that the corresponding value 
of B must be below critical. Combining both results gives 0.124568 < B^ < 0.124569, as 



announced in (3.2), and (1.2). 



Finally, dus B — Be is increased, the region over which the s ^ T^/^ scaling law holds 
becomes smaller and smaller, being overtaken by the s ^ law at high T joining the 
s ^ T law at low T, and ultimately will have disappeared altogether by the time B^ = 0.2 
is reached (this result is not shown in the figures). 




Figure 5: The left plot shows the crossover of s for low T. At moderately low temperatures 
s scales as T^/^ (1 ower left corner of the plot), while at ultra-low temperatures s scales as T 
for B > Be (curves a, b, c, d, e, f), and tends to a non-zero constant for B < B^ (curve g). 
The right plot shows the crossover for s from the moderately low temperature T^^^ scaling 
to the high temperature behavior. The dots represent numerical data points, while the 
solid interpolating lines are included to guide the eye. 
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3.5 Low B region 

Decreasing the magnetic field further, namely B < we find a completely different low 
temperature behavior: the entropy now stabilizes to a nonzero value at T = 0. We plot 
this limiting value in Figure [6j In pj] we speculated that an infinitesimally small magnetic 
field might be sufficient to remove the ground state entropy of the B = solution. We 
now see that this speculation is incorrect - only magnetic fields B > B^ accomplish this. 
The extremal entropy turns on continuously below the critical magnetic field. A curve that 
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Figure 6: Extremal entropy density: blue points represent results for the entropy density 
at very low temperatures in the region B < B^. The red curve is a plot of the function 

i75^^" (f ) 

passes through all the points to better than .5% accuracy is given by 



4^6 (1 + S2)7/4 



(3.4) 



The prefactor of 1/(4^/6) is fixed by the value for the purely electrically charged black brane, 
which is of course known analytically. It is remarkable that this function fits the data so 
well, but we are unable to say whether it represents an exact analytical result. 



To excellent accuracy, and as exhibited by (3.4), we find 



5- (3.5) 

near the critical point. If we increase the temperature at fixed magnetic field in this regime, 
we again find a region obeying the same s ^ T^/^ scaling as described in section 3.3. 
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3.6 Scaling region near the critical point 

The fact that the entropy density at the critical magnetic field has a simple power law 
dependence on temperature is indicative of a quantum critical point. Recall that since we are 
working in terms of dimensionless quantities as defined by the scalings in the asymptotically 
AdSs region, a priori the dimensionless entropy density s is allowed to be an arbitrary 
function of the dimensionless temperature T. This interpretation can be sharpened further 
by writing down a scaling form for the entropy in the vicinity of the critical point 

«^Ti/(^) (3.6) 

where, according to our results, the scaling function f{x) has asymptotic behavior 



^ < 




X — oc 

x^O (3.7) 

X ^ oo 



for some constants Ci^2,3- In general, the entropy is a function of two dimensionless combi- 
nations of 5, p, and T, whereas near the critical point the claim is that it can be written as 
a function of only one variable. All of this is of course standard from the general theory of 
classical and quantum critical phenomena. 

It is further useful to recall that for a quantum critical theory in d spatial dimensions, 
with dynamical critical exponent z, and with a relevant coupling m of scale dimension A, 
the entropy density will take the scaling form 



. = T-'^f ( (3.* 



Comparing to (3.7) we read off = 1, z = 3, and A = 2. 

We have verified that the scaling form ( |3.7[ ) conforms to our numerical results near the 
critical point B = Be and T = 0. From the numerics we can reconstruct the form of the 
scaling function, which is shown in Figure [7| We can use this to determine the constants 



appearing in (3.7) to be. 



ci?^ 0.172 C2?^ 0.110 0.045 (3.9) 

As noted in the Introduction, the standard approach to modelling a magnetically tuned 
quantum critical point of the type that we are seeing is based on the Hertz/Millis theory. 



We consider the action in (1.7) with d = 1. To the extent that this action captures the 
low energy degrees of freedom of our theory, it is natural to compute its finite temperature 
entropy density and compare to our results. This cannot be entirely correct for a number 



16 



fix) 



0.4- 




* 0.3 - 








0.1 - 

1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 



rnfHM{2,2x) 




X 



Figure 7: Plots of the scaling function /(x) controlling the thermodynamics within the scaling 
region. On the right plot we restrict to x > 0, corresponding to B > B^^ and compare the 
scaling function from gravity with a free field computation based on the Hertz/Millis theory; 
the latter is displayed as a red line. The factors of 1/10 and 2.2 in the latter are chosen to 
give a good match, and are expected since the the gravitational result contains a factor of 
Newton's constant and the normalizations of the magnetic fields need not agree. 

of reasons, not least that there is no way of explaining the ground state entropy density in 
this framework, but the comparison is instructive nonetheless. In the free field limit we just 



= B -Br 



need the dispersion relation implied by (1.7), which is 

uj=\k\{e + m^) 
The partition function is 

lnZ = -- r dklnil - e-^\'\^''+^'> 

271 J-oo 

from which we can extract the entropy density as 



(3.10) 



(3.11) 



(3.12) 



We can write the result in a scaling form as 



B - Be 

7^2/3 



(3.13) 
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and then compare with the scahng function coming from gravity. The comparison is shown 
in the right panel of Figure (rj We only compare in the region B > since it is only in this 
region that the Hertz/Millis action applies. Nothing fixes the normalization of the magnetic 
field appearing in ( |1.7D compared to that in gravity, and so we have allowed ourselves to 
adjust this relationship to achieve a good fit. Similarly, we have introduced a parameter 
to adjust the overall normalizations. The functions match surprisingly well, although it is 
unclear how much significance to attach to this. 

3.7 The quantum critical point 

Initial data yielding our closest approach to the quantum critical point is given by 

h = 1.7320507 

q = 5.08 X 10-^ (3.14) 
This choice of (6, along with the gauge choice C(r+) = 3.464063 leads to the following 



values for the asymptotic parameters of (2.8) 



V = 5.23811 X 10"^^ w = 1.60044 x 10"^ cq = 859.258 
63 = 5.20520 X 10^° ps = -2.08236 x 10^^ (3.15) 

from which we can compute the temperature and entropy density as 

f = 9.59602 X 10"^ 

s = 2.08699 X 10-^ (3.16) 



From (3.14) we see that there is a large hierarchy between the size of the electric and magnetic 
charges at the horizon in our chosen coordinates. Nevertheless, measured at infinity in the 
rest frame the ratio of these quantities is of order unity, B^ ^ 0.124569. This happens because 
there are large rescalings that occur in the process of integrating out from the horizon to the 
asymptotic regions, as illustrated by the magnitude of the parameters in (3.15). Using ( |2.9D 



these can convert a large hierarchy into one of order unity. 

It is instructive to examine the metric functions for this near critical solution. In Figure [8] 
we provide plots of these functions in the near (and not so near) horizon region. Near 
the horizon, one can think of these metric functions as representing a deformation of the 
configuration 

U{r) = 12(r-l)2 B = V3 

C{r) = 2A/3(r-l) V = W = E = P = (3.17) 

This represents AdSs x in "boosted coordinates" with magnetic fiux. We expect that. 



as we get nearer and nearer to the critical point, the metric functions will approach (3.17) 
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Figure 8: Plots of the metric functions for a solution very near the quantum critical point. 
Multiplicative factors have been included to allow the functions to appear on the same plot. 
The form of the functions near the horizon show that the metric is behaving as AdSa x 
plus small corrections. These small corrections are amplified in passing to the asymptotic 
AdSs regime. 



to greater and greater accuracy. On the other hand, it is crucial that there always be some 
deviation away from (3.17), as this solution is purely magnetic with vanishing electric charge 
density. Presumably what happens is that there is a fine-tuned limit in which all of the 
charge is carried by the bulk supergravity fields outside the near horizon region; this is 
possible by virtue of the Chern-Simons term. It would of course be extremely useful if this 
limiting solution could be found analytically. 



We close this subsection with an issue which has not yet been conclusively resolved by 
our numerical studies. For large 5, all fiows towards lower T end in the purely magnetic 
fixed point solution of [10], whose near-horizon geometry is AdS'^ x R^. More precisely, it 
can be shown numerically that all flows for B > end up at the purely magnetic fixed 
point. Where do the fiows for B < B^ end? Numerically, we have established that fiows for 
B less than but close to Be come near the purely magnetic fixed point, but at temperatures 
so low (on the order of T ^ 10~^^) that further numerical analysis can be pursued only at 
the cost of calculation times and memory which exceed those of our computers. Thus, the 
logical possibility that the fiows with B < Be end up on the critical curve, before the purely 
magnetic fixed point has been attained, cannot be excluded at this time. If this situation is 
in fact the one that occurs, then it would have a surprising similarity to the dynamics of the 
k = 1 theory, derived in [TT] . 
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3.8 Metamagnetic quantum criticality in Sr3Ru207 



We close this section with a few words about a real system that parallels ours in a number 
of ways. Quantum criticality and novel phases in Sr3Ru207 have been the subject of much 
experimental and theoretical interest in the past few years; e.g. [151 Ell- Sr3Ru207 is a 
layered structure, which for a large magnetic field perpendicular to the layers exhibits a 
line of first order metamagnetic phase transitions at finite temperature, ending at a finite 
temperature critical point. By including a component of magnetic field in the plane of the 
layers, the critical point can be brought to zero temperature. As in our case, the transition 
occurs at finite magnetic field, and involves no change of symmetry. Away from the critical 
point the system behaves as a Fermi liquid, with entropy density linear in temperature. As 
the critical point is approached, the linear term diverges as s/T ^ 1/{B — B^)^ just as we 
found (in this case is approximately 8 Tesla.) In recent work, the complete "entropy 
landscape" of Sr3Ru207 at finite temperature and magnetic field has been mapped out [15]. 
In very pure samples, as one tries to sit right on top of the critical point one finds instead 
that a new phase emerges, which is believed to be a spatially anisotropic nematic phase 
f33l [34] . This has been described as nature's solution to the problem of avoiding a finite 
entropy density at zero temperature. The parallels with our system are evident (though an 
obvious diflFerence is that our critical theory eflFectively sees only a single spatial direction, 
compared to the two in-plane directions in Sr3Ru207), and lead us to speculate whether 
in our case the extremal black hole phase is unstable and gives way to a new phase with 
reduced symmetry, as in [35j. 



4 Discussion 

We have found a holographic description of a quantum critical point, reached by tuning a 
magnetic field at finite density, that nicely resembles examples seen in the real world. In 
particular, approached from the high 5-field side, we have Fermi liquid behavior; as the 
magnetic field is lowered the specific heat coefficient diverges, and we enter a regime of 
non- Fermi liquid behavior with nontrivial scaling properties. 

The most exotic property as compared to known physical examples is the existence of 
a ground state entropy on the low field side of the transition. The presence of a non- 
vanishing entropy at T = for 5 < 5c would again seem to contradict the third "law" 
of thermodynamics [H], as it did in the absence of magnetic fields. Actually, our results 
show that the lifting of the ground state entropy as a function of an external magnetic field 
is a subtle and dynamical issue. Also, we expect the ground state entropy to be lifted by 
instabilities driving the system towards inhomogeneity, or by turning on further fields besides 
an external magnetic field. 

It is worth emphasizing again the universality of our results: they apply to all supersym- 
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metric AdSs examples arising from IIB or M-theory, since all such theories admit a consistent 
truncation to the Einstein-Maxwell-Chern-Simons theory used here (UEJIS]. We did not have 
to introduce any model building devices in the way of probe branes or scalar fields. 

There are many open questions and avenues for further investigation. Many of our 
numerical results cry out for an analytical derivation. In particular, one would expect to be 
able to derive the value z = 3 of the dynamical exponent, and the dimension A = 2 of the 
relevant operator governing the critical theory. It may similarly be possible to understand 
these results microscopically on the gauge theory side. 

All of our results were presented for the supersymmetric value of the Chern-Simons 
coupling, k = 2/^/3, but it would be useful to consider other values as well. Our expectation 
is that as k is increased the critical point will retain its character but with Be moving to a 
smaller value. An interesting question would then be whether we reach 5c = at finite fc, 
for if so there would no longer be an extremal black hole phase. 

In this work we have only considered the thermodynamics, but the existence of the 
critical theory implies scaling behavior of correlation functions with respect to frequency and 
momentum. Computing these would be valuable in pinning down the precise connection to 
the Hertz/Millis theory. 

Finally, to more closely model real systems it would be very interesting to construct a 
version of our system giving rise to critical behavior in two or three spatial dimensions. 
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